Nonequilibrium relaxation of the two-dimensional Ising model: Series-expansion and 

Monte Carlo studies 
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We study the critical relaxation of the two-dimensional 
Ising model from a fully ordered configuration by series ex- 
pansion in time t and by Monte Carlo simulation. Both the 
magnetization (m) and energy series are obtained up to 12- 
th order. An accurate estimate from series analysis for the 
dynamical critical exponent z is difficult but compatible with 
2.2. We also use Monte Carlo simulation to determine an 
effective exponent, z^sit) = —^d\nt/d\nm, directly from a 
ratio of three-spin correlation to m. Extrapolation to t ^ oo 
leads to an estimate z = 2.169 ± 0.003. 

PACS number(s): 05.50.-hq, 05.70.Jk, 02.70.Lq 



I. INTRODUCTION 

The pure relaxational dynamics of the kinetic Ising 
model with no conserved fields, which is designated as 
model A in the Hohenberg-Halperin review [Q , has been 
studied extensively by various approaches. Unhke some 
of the other models in which the dynamical critical ex- 
ponent z can be related to the static exponents, it seems 
that z of model A is independent of the static exponents 
(however, see Rcf. Q). In the past twenty years, the 
numerical estimates for the dynamical critical exponent 
z scattered a lot, but recent studies seem to indicate a 
convergence of estimated values. Our studies contribute 
further to this trend. 

We review briefly some of the previous work on the 
computation of the dynamical critical exponents, concen- 
trating mostly on the two-dimensional Ising model. The 
conventional theory predicts z — 2 — rj |^ , where rj is the 
critical exponent in the two-point correlation function, 
G{r) (X For the two-dimensional Ising model, 

this gives z = 1.75. It is known that this is only a lower 
bound H- It is very interesting to note that series expan- 
sions p^TlOl] gave one of the earliest quantitative estimates 
of z. Dammann and Reger |l^ have the longest high- 
temperature series (20 terms) for the relaxation times so 
far, obtaining z — 2.183 ± 0.005. However, re-analysis of 
the series by Adler gives z = 2.165±0.015. There are 
two types of field-theoretic renormalizatio n group analy- 
sis: the e-expansion near dimension d = 4 |12| , p^ and an 
interface model near d = 1 |l^. It is not clear how reli- 
able when it is interpolated to d = 2. Real-space renor- 



malization group of various schemes has been proposed in 
the early eighties ||l^-|l^, but it appears that there are 
controversies as whether some of the schemes are well- 
defined. The results are not of high accuracy compared 
to other methods. Dynamic Monte Carlo renormaliza- 
tion group p9|-|2^ is a generalization of the equilibrium 
Monte Carlo renormalization group method [2^. The 
latest work |2^ gives z — 2.13 ± 0.01 in two dimensions. 
Equilibrium Monte Carlo method is one of the standard 
methods to estimate z [|2^-|2^. However, long simula- 
tions {t ^ L^) are needed for sufficient statistical ac- 
curacy of the time-displaced correlation functions. The 
analysis is quite difficult due to unknown nature of the 
correlation functions. Nonequilibrium relaxation [p9|-p4|, 
starting from a completely ordered state at Tc, has nice 
features. The analysis of data is more or less straightfor- 
ward. The lattice can be made very large, so that finite- 
size effect can be ignored (for t ^ L^). The catch here 
is that correction to scaling due to finite t is large. Re- 
cently, the idea of damage spreading 135^391 has also been 
employed. Methods based on statistical errors in equi- 
librium Monte Carlo simulation jEo| , finite-size scaling 
of nonequilibrium relaxation pl| , [42| , and finite-size scal- 
ing of the eigenvalues of the stochastic matrix are 
used to compute the exponent. A recent calculation with 
a variance-reducing Monte Carlo algorithm for the lead- 
ing eigenvalues gives prediction Q z — 2.1665 ± 0.0012. 
This appears to be the most precise value reported in the 
literature. 

The high-temperature series expansions for the relax- 
ation times are often used in the study of Ising dynamics. 
In this paper, we present a new series which directly cor- 
responds to the magnetization (or energy) relaxation at 
the critical temperature. Our series expansion method 
appears to be the only work which uses time t as an 
expansion parameter. The generation of these series is 
discussed in Sections |l| and HI. Dynamical scaling men- 
tioned in Section |^ forms the basis of the analysis, and 
the results are analyzed in Section We feel that the 
series are still too short to capture the dynamics at the 
scaling regime. We also report results of an extensive 
Monte Carlo simulation for the magnetization relaxation. 
We find that it is advantageous to compute an effective 
dynamical critical exponent directly with the help of the 
governing master equation (or the rate equation). The 
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simulation and analysis of Monte Carlo data are pre- 
sented in Section VI, We summarize and conclude in 



Section VII 



II. SERIES EXPANSION METHOD 

In this section, we introduce the relevant notations, 
and outline our method of series expansion in time vari- 
able t. The formulation of single-spin dynamics has al- 
ready been worked out by Glauber and by Yahata 
and Suzuki |^ long time ago. To our knowledge, all 
the previous series studies for Ising dynamics |5|-|lC|| are 
based on high-temperature expansions of some correla- 
tion times. As we will see, expansion in t is simple in 
structure, and it offers at least a useful alternative for 
the study of Ising relaxation dynamics. 

We consider the standard Ising model on a square lat- 
tice [Eg! with the energy of a configuration a given by 



E(a) 



-J ^ a^aj, 



(1) 



where the spin variables ai take ±1, J is the coupling 
constant, and the summation runs over all nearest neigh- 
bor pairs. The thermal equilibrium value of an observ- 
able /(cr) at temperature T is computed according to the 
Boltzmann distribution. 



Y.aexp{-E{a)/kBT) 



^/(a)Peq(a). 



(2) 



The equilibrium statistical-mechanical model defined 
above has no intrinsic dynamics. A local stochastic dy- 
namics can be given and realized in Monte Carlo simula- 
tions |4^ . The dynamics is far from unique; in particular, 
cluster dynamics differs vastly from the local ones. 

A sequence of Monte Carlo updates can be viewed as a 
discrete Markov process. The evolution of the probability 
distribution is given by 



P(a,fc + l) = ^P(a',fc)W^(aV), 



(3) 



where is a transition matrix satisfying the stationary 
condition with respect to the equilibrium distribution, 
i.e., Poq ~ PcqW. A continuous time description is more 
convenient for analytic treatment. This can be obtained 
by fixing t — k/N, and letting 6t = 1/N — 0, where N = 
is the number of spins in the system. The resulting 
differential equation is given by 



dP{a,t) 
dt 



= TP{a,t), 



(4) 



where F is a linear operator acting on the vector P((T, t), 
which can be viewed as a vector of dimension 2^, indexed 



by cr. If we use the single-spin-flip Glauber dynamics [Q 
we can write 



N N 

j=i 3=1 



where 



1 — (Jj tanh^/^ (Tfc 

nil of j 



K = 



knT' 



and Fj is a flip operator such that 



= P(...,-a„...). 



(5) 



(6) 



(7) 



The flip rate Wj {aj ) for site j depends on the spin value 
at the site j as well as the values of its nearest neighbor 
spins CTfc. 

The full probability distribution clearly contains all the 
dynamic properties of the system. Unfortunately its high 
dimensionality is difficult to handle. It can be shown 
from the master equation, Eq. (^), that any function of 
the state a (without explicit t dependence) obeys the 
equation 



d{f)t 
dt 



where 



N 



and the average of / at time t is defined by 



(8) 



(9) 



(10) 



Note that the time dependence of (/) is only due to 
P{a,t). For the series expansion of this work, it is suf- 
ficient to look at a special class of functions of the form 
= Wj^A ^3^ where A is a set of sites. In such a case 
we have 



dt 



(11) 



With this set of equations, we can compute the n-th 
derivative of the average magnetization (ao)t. A formal 
solution to Eq. (H) is 



n=0 



(12) 



This equation or equivalently the rate equation, Eq. ([ll]), 
forms the basis of our series expansion in time t. 

A few words on high-temperature expansions are in 
order here. They are typically done by integrating out 
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the time dependence — the nonhnear relaxation time can 
be defined as 



/o 



(13) 

The equihbrium correlation time (linear relaxation time) 
can be expressed as 



(m(0)m(t))cq V- 



\"l(0)2)cq ^" ^ 



(14) 



where x — N{m'^)cq is the reduced static susceptibility. 
The average is with respect to the equilibrium distribu- 
tion, Pcq{a). A suitable expansion in small parameter 
J/kBT can be made by writing C = Co + AC. 

It is clear that we can also perform the Kawasaki dy- 
namics with a corresponding rate. Of course, since the 
magnetization is conserved, only energy and higher order 
correlations can relax. 

A very convenient form for the Glauber transition rate, 
Eq. (^, on a two-dimensional square lattice is 

+yao{ai(T2a3 + a2(T3a4 + o'3<740'i + (y4(7i(T2) , (15) 



X = -- tanh 2K -- tanh iK, 
4 8 



y = +- tanh 2K - - tanh iK, 



(16) 



(17) 



where the site is the center site, and sites 1, 2, 3, and 4 
are the nearest neighbors of the center site. At the critical 
temperature, taiihKc = v^— 1, we have X = -5\/2/24 
and y = \/2/24. 



III. COMPUTER IMPLEMENTATION AND 
RESULTS 

A series expansion in t amounts to finding the deriva- 
tives evaluated at t = 0: 



E 

ri=0 



i! dt" 



t=0 



(18) 



The derivatives are computed using Eq. (|11|) recursively. 
A general function is coded in C programming language 
to find the right-hand side of Eq. ( pi] ) when the configu- 
ration a^, or the set A, is given. The set A is represented 



as a list of coordinates constructed in an ordered man- 
ner. By specializing the fiip rate as given by Eq. (|l^), 
and considering each site in A in turn, the configurations 
on the right-hand side of the rate equation are generated 
in three ways: (1) the same configuration as A, which 
contributes a factor (coefficient of a term) of —1; (2) a 
set of configurations generated by introducing a pair of 
nearest neighbor sites in four possible directions, with one 
of the sites being the site in A under consideration, and 
making use of the fact af = 1. We notice that the site 
in A under consideration always gets annihilated. Each 
resulting configuration contributes a factor of —x; and 
(3) same as in (2) but two more sites which are also the 
nearest neighbors of the site in A under consideration are 
introduced. This two extra sites form a line perpendicu- 
lar to the line joined by the first pair of neighbor sites in 
(2). Each of this configuration has a factor of —y. It is 
instructive to write down the first rate equation, taking 
into account of the lattice symmetry (e.g. (ai) — (ctq), 
for all i): 



d{(To) 
dt 



= -(1 + 4a;)(o-o) - 4y(cricr2CT3) 



(19) 



The core of the computer implementation for series 
expansion | |49[ | is a symbolic representation of the rate 
equations. Each rate equation is represented by a node 
together with a list of pointers to other nodes. Each 
node represents a function (cr^), and is characterized by 
the set of spins A. The node contains pointers to the 
derivatives of this node obtained so far, and pointers to 
the "children" of this node and their associated coeffi- 
cients, which form a symbolic representation of the rate 
equations. The derivatives are represented as polynomi- 
als in y. Since each node is linked to other nodes, the 
computation of the n-th derivative can be thought of as 
expanding a tree (with arbitrary number of branches) of 
depth n. 

The traversal or expansion of the tree can be done in 
a depth-first fashion or a breadth-first fashion. Each has 
a different computational complexity. A simple depth- 
first traversal requires only a small amount of memory of 
order n. However, the time complexity is at least expo- 
nential, 6", with a large base b. A breadth-first algorithm 
consumes memory exponentially, even after the number 
of the rate equations has been reduced by taking the 
symmetry of the problem into account. The idea of dy- 
namic programming can be incorporated in the breadth- 
first expansion where the intermediate results are stored 
and referred. To achieve the best performance, a hybrid 
of strategies is used to reduce the computational com- 
plexity: 

• Each configuration (pattern) is transformed into 
its canonical representation, since all configurations 
related by lattice symmetry are considered as the 
same configuration. 

• We use breadth-first expansion to avoid repeated 
computations involving the same configuration. If 
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a configuration has already appeared in earlier ex- 
pansion, a pointer reference is made to the old con- 
figuration. Each configuration is stored in memory 
only once. However, storing of all the distinct con- 
figurations leads to a very fast growth in memory 
consumption. 

• The last few generations in the tree expansion use 
a simple depth-first traversal to curb the problem 
of memory explosion. 

• Parallel computation proves to be useful. The 
longest series is obtained by a cluster of 16 Pen- 
tium Pro PCs with high speed network connection 
(known as Beowulf). 

The program is controlled by two parameters D and 
C. 13 is the depth of breadth- first expansion of the tree. 
When depth D is reached, we no longer want to continue 
the normal expansion in order to conserve memory. In- 
stead, we consider each leave node afresh as the root of 
a new tree. The derivatives up to (n — D)th order are 
computed for this leave node. The expansion of the leave 
nodes are done in serial, so that the memory resource can 
be reused. The parameter C controls the number of last 
C generations which should be computed with a simple 
depth-first expansion algorithm. It is a simple recursive 
counting algorithm, which uses very little memory, and 
can run fast if the depth C is not very large. In this 
algorithm the lattice symmetry is not treated. The best 
choice of parameters is 13 = 6 and C = 2 on a DEC 
AlphaStation 250/266. The computer time and mem- 
ory usage are presented in Table m. As we can see from 
the table, each new order requires more than a factor of 
ten CPU time and about the same factor for memory if 
memory is not reused. This is the case until the order 
D + C + 1, where no fresh leave- node expansion is made. 
There is a big jump (a factor of 60) in CPU time from 9- 
th order to 10-th order, but with a much smaller increase 
in memory usage. This is due to the change of expan- 
sion strategy. Finally the longest 12-th order series is 
obtained by parallel computation on a 16-node Pentium 
Pro 200 MHz cluster in 12 days. The number of distinct 
nodes generated to order n is roughly ^^11". To 12-th 
order, we have examined about 10^*^ distinct nodes. The 
series data are listed in Tabic ||. 



IV. DYNAMICAL SCALING 

The traditional method of determining the dynami- 
cal critical exponent z is to consider the time-displaced 
equilibrium correlation functions. However, one can al- 
ternatively look at the relaxation towards thermal equi- 
libration. The basic assumption is the algebraic decay of 
the magnetization at Tc, 



ct 



-13/i^z 



This scaling law can be obtained intuitively as follows. 
Since the relaxation time and the correlation length are 
related through t cx by definition, after time t, the 
equilibrated region is of size i^/^. Each of such region 
is independent of the others, so the system behaves as 
a finite system of linear length ^ (x t^^^. According to 
finite-size scaling [Q, the magnetization is of order ^"^Z" 
on a finite system of length ^. Each region should have 
the same sign for the magnetization since we started the 
system with all spins pointing in the same direction. The 
total magnetization is equal to that of a correlated region, 
giving m (x t~^^^^ . 

The same relation can be derived from a more general 
scaling assumption pit]. 



T-T, 



(21) 



By requiring that m{t, e) is still finite as the scalin g ar - 
gument te'^^ and e ^ with fixed t, we get Eq. ([2(]). 

Equation ( |20| ) is only true asymptotically for large t. It 
seems that there is no theory concerning leading correc- 
tion to the scaling. As a working hypothesis, we assume 
that 



m 



ct-'^/'^'il + bt-^). 



(22) 



The Monte Carlo simulation results as well as current 
series analysis seem to support this with A near 1. Other 
pos sibility might he z — 2 with logarithmic correction 



V. ANALYSIS OF SERIES 

A general method for extending the range of conver- 
gence of a series is the Pade analysis |53 5^ where a series 
is approximated by a ratio of two polynomials. We first 
look at the poles and zeros of the Pade approximants in 
variable s — t/{t + 1) for m. Since t varies in the range 
of [0, oo), it is easier to look at s, which maps the inter- 
val [0, oo) to [0, 1). There are clusters of zeros and poles 
in the s-interval (1,2) which corresponds to negative t. 
But interval [0, 1) is clear of singularities, which gives 
us hope for analytic continuation to the whole interval 
[0, 1). If we assume the asymptotic behavior m oc i^", 
then dlnm/dt = —a/t « — a(l — s) for large t or s — > 1. 
This means that the Pade approximant should give a zero 
around s = 1. We do observe zeros near 1. But it is typ- 
ically a pair of zeros off the real axis together with a pole 
at the real axis near 1, or sometimes, only a pair of real 
zeros. These complications make a quantitative analysis 
difhcult. 

Since wc know the exact singular point (corresponding 
to t — oo), we use the biased estimates by considering 
the function 



t 



(20) 



Fit) = 



dhim 
dint 



A 

vz 



(23) 
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An effective exponent Zcs{t) is defined by Zcff(i) = 
-P/{„Fit)) = -l/{8Fit)). 

Again we prefer to use tlie variable s to bring the in- 
finity to a finite value 1. Due to an invariance theorem 
, the diagonal Pade approximants in s and t are equal 
exactly. For off-diagonal Pade approximants, s is more 
useful since the approximants do not diverge to infinity. 

We use methods similar to that of Dickman et al. pq ] 
and Adler The general idea is to transform the 

function 171(1) into other functions which one hopes to be 
better behaved than the original function. In particular, 
we require that ast—^ 00, the function approaches a con- 
stant related to the exponent z. The first transformation 
is the Eq. (pi) . A second family of transformations is 



Gpit) 



d\nJ^m{t')Pdt' 



dint 



8z' 



(24) 



where p is a real positive number. One can show that the 
two functions are related by 



The last transform is 



(25) 



(26) 



where A is an adjustable parameter, and F can also be 
replaced by Gp. If the leading correction to the constant 
part is of the form , the transformation will eliminate 
this correction term. 

The transformation of the independent variable t to 
other variable is important to improve the convergence 
of the Pade approximants. We found that it is useful to 
consider a generalization of the Euler transform. 



1 



1 



(1 + *)^ 



(27) 



The parameter A is adjusted in such a way to get best 
convergence among the approximants. Since for i — > 00 
or u 1, a Pade approximant near m = 1 is an analytic 
function in u, which implies that the leading correlation 
is of the form t~'^. Note that A = 1 corresponds to the 
Euler transformation [u = s when A = 1). 

One of the fundamental difficulties of the transforma- 
tion method is that one does not know a priori that 
a certain transformation is better than others. Worst 
still, we can easily get misleading apparent convergence 
among different approximants. Thus, we need to be very 
careful in interpreting our data. Specifically, we found 
that Eq. (|2^) gives less satisfactory result than that of 
Eq. (|2^) , where the independent variable t is transformed 
into u according to Eq. (^7|) . Figure |l| is a plot of all the 
Pade approximants of order [N, D], with N > A, D > A, 
and TV -I- Z) < 12, as a function of the parameter A, for 
Gi{t = 00). Good convergence is obtained at A = 1.217 
with z « 2.170. The estimates z vary only slightly with 



p, at about 0.005 as p varies from 0.5 to 2. Using F{t) 
of Eq. (H), the optimal value is A = 1.4 with z « 2.26. 
Using the function H does not seem to improve the con- 
vergence. Even though the value 2.170 seems to be a 
very good result, we are unsure of its significance since 
there are large deviations of the Pade approximation to 
the function F{t) for 1/t < 0.2 from the Monte Carlo 
result of Fig. ||. 

An objective error estimate is difficult to give. Esti- 
mates from the standard deviation of the approximants 
tend to give a very small error but incompatible among 
different methods of analysis. Different Pade approx- 
imants are definitely not independent; we found that 
[N, D] Pade is almost equal to [D, N] Pade to a high 
precision. A conservative error we quote from the series 
analysis is 0.1. 

Analysis of the energy series is carried out similarly 
with m replaced by (crocri) — ^2/2, where the constant 
V2/2 is the equilibrium value. The large t asymptotic 
behavior is t~^/^ Both F and G functions give com- 
parable results, better convergence is obtained for A > 1. 
The value for z is about 2.2, but good crossing of the 
approximants are not observed. We feel better analysis 
method or longer series is needed. 



VI. MONTE CARLO SIMULATION 

Our motivation for a Monte Carlo calculation was to 
check the series result. It turns out that the data are 
sufficiently accurate to be discussed in their own right. 
Such an improved accuracy is achieved by using Eq. (|lS|), 
which permits a direct evaluation of the effective expo- 
nent Zcs{t). 

We compute the magnetization m = (cro): energy 
per bond (croo'i), and the three-spin correlation 7713 — 
{(Jia2(T3) where the three spins are the nearest neighbors 
of a center site having one of the neighbor missing in the 
product. With these quantities, the logarithmic deriva- 
tive, Eq. (|23|), can be computed exactly without resort- 
ing to finite differences. From Eq. ([l9| ) we can write (at 

T^Tr) 



„ , , d In m 

= 7h^ — ^ 



V2 



6 v TO 



8Zeff(0 ■ 

(28) 



The above equation also defines the effective exponent 
Zcs{t) which should approach the true exponent z as t — > 
00. 

The estimates for the effective exponent based on the 
ratio of one spin to three-spin correlation, Eq. (p8|), have 
smaller statistical errors in comparison to a finite differ- 
ence scheme based on m{t) and m{t + l). Error propaga- 
tion analysis shows that the latter has an error 5 times 
larger. Both methods suffer from the same problem that 
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error 5z oc t. Thus, working with very large t does not 
necessarily lead to any advantage. 

In order to use Eq. (^Sj), we need exactly the same flip 
rate as in the analytic calculations, namely the Glauber 
rate, Eq. (^. The continuous time dynamics corresponds 
to a random selection of a site in each step. Sequential 
or checker-board updating cannot directly be compared 
with the analytic results. However, it is believed that 
the dynamical critical exponent z does not depend on 
the details of the dynamics. 

We note that a Monte Carlo simulation is precisely 
described by a discrete Markov process while the series 
expansion is based on the continuous master equation. 
However, the approach to the continuous limit should 
be very fast since it is controlled by the system size — 
the discreteness in time is We have used a sys- 

tem of 10^ X 10^, which is sufficiently large. Apart from 
the above consideration, we also checked finite-size effect. 
Clearly, as t > L^, finite-size effect begins to show up. 
We start the system with all spins up, m(0) = 1, and 
follow the system to t = 99. For t < 100, we did not 
find any systematic finite-size effect for L > 10'^. So the 
finite-size effect at L = 10^ and t < 100 can be safely 
ignored. 

Figure || shows the Monte Carlo result for the effective 
exponent as a function of 1/t. The quantities m, 7713, 
and {(Joai) are averaged over 1868 runs, each with a sys- 
tem of 10^ spins. The total amount of spin updating is 
comparable to the longest runs reported in the literature. 
Based on a least-squares fit from t = 30 to 99, we obtain 



the high-temperature series. We have also studied the re- 
laxation process with Monte Carlo simulation. The ratio 
of three-spin to magnetization is used to give a numeri- 
cal estimate of the logarithmic derivatives directly. This 
method gives a more accurate estimate for the dynamical 
critical exponent. 
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z = 2.169 ±0.003. 



(29) 



The error is obtained from the standard deviation of few 
groups of independent runs. An error estimate based on 
the residues in the linear least-squares fit is only half of 
the above value, which is understandable since the points 
in Fig. ^ are not statistically independent. 

In Fig. ^, we also plot a series result for the Fit), 
obtained from the [6,6] Fade of Gi{u) and Eq. (25). 



Substantial deviations are observed for 1/t < 0.2, even 
though in the 1/t ^ limit, both results arc almost the 
same. This casts some doubts on the reliability of the 
series analysis. We note that the t ^ 00 limit of the 
function F{t) is invariant against any transformation in 
t which maps i = cxd to cxd. Thus, the discrepancy might 
be eliminated by a suitable transformation in the Fade 
analysis. 



VII. CONCLUSION 



We have computed series for the relaxation of mag- 
netization and energy at the critical point. The same 
method can be used to obtain series at other tempera- 
tures or for other correlation functions. The analyses of 
the series are non-trivial. We may need much more terms 
before we can obtain result with accuracy comparable to 
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TABLE I. CPU time and memory usage for the series ex- 
pansion of relaxation of magnetization, measured on an Al- 
phaStation 250/266. 



n CPU time (sec) Memory (MB) 

~6 aiS 003 

7 1.8 0.27 

8 25 3 

9 358 34 

10 23600 51 

11 939000 70 
12=" 1.6 X 10^ 85 



^Actual computations are done on a 16-node Pentium Pro 200 
cluster. 



TABLE IL Series-expansion coefficients (n-th derivative) for a single spin {cTo)t and nearest neighbor spin correlation {aoai)t. 



n 


£i"(o-o>t 


rft" 1" 





1 


1 


1 


-1 + {2V2)/S 


-2 + (4^/2)/3 


2 


13/9 -%/2 


(56 - 39^/2)/9 


3 


(15- ll^/2)/27 


2(-249 + 175\/2)/27 


4 


-53/3 + 25/V2 


(1988 - 1399\/2)/54 


5 


(41175 - 29111\/2)/486 


(30834 - 21919\^)/486 


6 


(-66133 + 46680i^)/1458 


2(-142869 + 101087\/2)/243 


7 


(-125718825 + 88903747v^)/34992 


5(18191091 - 1286740iy2)/17496 


8 


17(92513582 - 6541830iy2)/34992 


(2190719830 - 1548846809 V2)/69984 


9 


(-429437553903 + 303660675715V2)/1259712 


(-289028693217 + 204371 192813V2)/314928 


10 


(4931635327666 - 3487215692619a/2)/3779136 


(43146864055759 - 30509318092215\/2)/3779136 


11 


(1821425391381531 - 1287938652305897V2)/181398528 


(-957792089655213 + 677259915390707 V2)/10077696 


12 


7(-10761633667757321-|-7609621330268025a/2)/272097792 


(425962164223774298-301200006005168631^2) /1088391 168 



FIG. 1. Pade estimates of the dynamical critical exponent 
z using Gi{t — (x), plotted as a function of A, the transfor- 
mation parameter. On this scale, the Pade approximant of 
order [N,D] is indistinguishable from [D,N]. 



FIG. 2. Effective exponent 2cff(t) plotted against inverse 
time 1/ t. The circles are Monte Carlo estimates based on 
Eq. (pq); the continuous curve is obtained from the [6,6] 
Pade approximant of Gi in variable u, transformed back to 
F through Eq. (pHI). 
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Fig. 2, Wang & Gan 



Fig.l, Wang & Gan 




